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Abstract 

Chirped Bessel waves are introduced as stable (non-diffracting) solutions of the paraxial wave 
equation in optical antiguides with a power-law radial variation in their index of refraction. 
Through numerical simulations, we investigate the propagation of apodized (finite-energy) ver¬ 
sions of such waves, with or without vorticity, in antiguides with practical parameters. The new 
waves exhibit a remarkable resistance against the defocusing effect of the unstable index potentials, 
outperforming standard Gaussians with the same full width at half maximum. The chirped prohle 
persists even under conditions of eccentric launching or antiguide bending and is also capable of 
self-healing like standard diffraction-free beams in free space. 


* ioannis.chreminos@mpl.mpg.de 


1 




I. INTRODUCTION 


Structured light waves have become particularly interesti^ in recent years due to their 
remarkable properties and extended range of applications. [l(] A big part of the effort in 
the held has been targeted at the so-called non-diffracting waves, namely waves that re¬ 
sist diffraction (in their realistic hnite-energy versions) and thus maintain their transverse 

n n 

intensity prohle over large distances. [21] Bessel beams [3| are perhaps the most familiar ex¬ 
ample with diverse applications in optical manipulatiom atom and nonlinear optics, jd] 


Other examples are Mathieu [^ and parabolic beams jo] which follow as separable solu¬ 
tions of the wave equation in elliptic or parabolic coordinates, respectively. Self-imaging 
beams that reproduce their prohle periodically with distance and beams that rotate their 
prohle around the axis of pr^agation also qualify as non-dihracting or, more generally, as 
propagation-invariant waves. [7j More recently, the family of non-dihracting beams expanded 
with the introduction of accelerati^ Airy beams js] and, subsequently, with the complete 
set of paraxial accelerating beams. |9| These waves have a non-dihracting prohle that shifts 
laterally according to a quadratic law with the propagation distance, thus writing a parabola 
in free space. The demonstrated and envisaged applications of accelereting beams have also 
been numerous. 0 Hybrids between the accelerating and non-accelerating families are also 
oossible, as for example the accelerating Bessel-like beams in the paraxial [ll| or nonparaxial 
12| regime. 

Non-dihracting waves owe their properties to their special structure of rays which is 
such that the resulting transverse interference pattern does not change with the propagated 
distance. A Bessel beam, for example, is produced Iw ray cones with a hxed opening angle 
and expanding circular bases on the input aperture, (stj With appropriate deformation of the 
cones, their apexes can be made to trace a curved path which results in an accelerating 
Bessel-like beam. [ m| Airy beams, on the other hand, have a very diherent structure since 
their rays create parabolic optical caustics. What is however common among non-dihracting 
waves is that their prohle at increasing propagation distances is produced by rays emitted 
from points on the input aperture that are increasingly far from their “center of mass” (The 
center of mass of a beam with amplitude u propagating along the ^ axis lies at the transverse 
position (J r\u\'^dxdy'j / (/ \u\‘^dxdy)). The on-axis far held of a Bessel beam, for example, 
is produced by rays emitted from the expanding circles on the input plane mentioned before. 
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This property of non-diffracting beams ensures that, if some of the inner rays are scattered 


or blocked, the outer rays will still produce the beam’s prohle at longer c 


is a remarkable property of non-diffracting beams called self-healing. 13|, [l^ A wave-optics 


istances, which 


description of this property is also possible. 1^ This should be contrasted with Gaussian 
beams where the on-axis field is produced by the single on-axis ray (except on the focal 
point). If this ray is blocked then the field along the entire z axis is distorted. Along the 
same lines, certain advanta ges have been highlighted in the propagation of non-diffracting 

le-Esl 


beams in turbulent media. 


Perusing the literature one sees that non-diffracting waves have been so far a concept 
limited in free space or homogeneous media, with medium inhomogeneities introduced only 
locally or randomly when self-healing after scattering or robustness in turbulence are to be 
studied, respectively. Going beyond these ideas, one may wish to explore the possibility of 
waves that propagate in a non-diffracting fashion in non-scattering inhomogeneous media 
that are inherently more diffractive than free space. That would be for example the case of 
optical antiguides (also called unstable optical ducts 19|), namely cylindrically symmetric 
inhomogeneous media whose index of refraction increases with the distance from the axis 
thus making rays bend outwards and waves diffract quickly. 

That non-diffracting waves of this kind should exist is already suggested by the Airy 
wavepacket paradigm. The Airy wavefunction is the unique non-spreading solution (exclud¬ 
ing the unbounded function Bi) of the one-dimensional (ID) Schrodinger equation in a field 
of constant force, By the analogy between Schrodinger equation and Helmholtz equation 
for paraxial optical waves, it is deduced that a ID optical Airy beam propagates without 
diffracting (and without accelerating!) in a medium whose index of refraction varies linearly 
in the transverse direction, e.g. n{x) = no -|- ax. This is obvious from the fact that the Airy 
function solves the differential equation u"{x) = xu{x). Again, the non-spreading nature 
(or stability) of this wave can be explained by its structure of parabolic rays enveloping a 


straight caustic. 


20l | Now imagine a cylindrical antiguide whose index increases linearly with 


the distance p from the axis. If we consider circularly symmetric waves, their corresponding 
rays are meridional and the underlying ray equation is identical to the ID case with the 
cartesian coordinate x replaced by the radial coordinate p. Hence, a circularly symmetric 
beam must exist (a 2D wave) that propagates without diffracting along the unstable linear 
duct and has in the p — z plane the same ray structure with the ID Airy wavepacket in a 
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linear potential. 

Beyond these simple arguments, in this paper we report a family of waves that propagate 
without diffracting in optical antiguides with a power-law variation of their refractive index. 
These waves can have zero or nonzero angular momentum while their radial amplitude 
satishes a Bocher-type differential equation. At sufficiently large radial positions the solution 
behaves asymptotically like a Bessel function with an argument of the order O , 

where u is the power with which the refractive index varies. For a certain value of the 
propagation constant that is equal to the wave number on the antiguide’s axis (core wave 
number), the wave amplitude is exactly proportional to a Bessel function with an argument 
that is proportional to Jn any case, this nonlinear dependence of the argument 

on the radius causes the wave amplitude to oscillate with an increasing spatial frequency 
similar to a chirped pulse in time. These properties prompt us to name the new waves chirped 
Bessel waves. Similar to standard Bessel beams in homogeneous space, chirped Bessel waves 
in antiguides carry ideally inhnite power which means that their realistic versions have to 
be subject to some kind of apodization of aperturing. 

Closing this introduction, we would like to mention that the question of trapping vortex 


waves in an optical antiguide was addressed recently. 


21| In this work simple vortex modes 


of the Laguerre-Gauss type were launched in curved antiguides with a Gaussian index prohle 
hoping to observe some trapping effect and guidance of the vortex along the antiguide. The 
conclusion was negative as the test waves had no special structure (for example chirp) in the 
radial direction. The chirped Bessel waves introduced here can indeed propagate in a stable 
(diffraction-resisting) fashion along a straight optical antiguide. Moreover, our simulations 
show that these waves remain robust when launched in a curved antiguide being however 
attracted along a curved trajectory toward the direction of increasing index. 

II. THEORY 

A. Wave analysis 


To keep things simple, we consider paraxial propagation of waves in an optical antiguide 
with a refractive index that increases radially as 

1/2 


n{p) = no 


1 + 2A 


P 


( 1 ) 
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where a is a characteristic radius and A > 0, z/ > 1 are real parameters. Notice that in 
(fTj) we have used a notation that is familiar within the held of graded-index optical hbers. 22| 
To stay in the paraxial regime and ignore polarization effects, the increase in the refractive 
index must act as a perturbation to the core index hq. In analogy to the terminology of 
weakly guiding waveguides and hbers, an antiguide of this form can be thought to be weakly 
antiguiding. With A 1, this will be true within distances of few a from the axis. We will 
return to this point later in this section. 

In the paraxial regime, any component of the electromagnetic held can be expressed as 
n(p, (p, where u is an envelope function that varies slowly along and fco = ‘^ttuq/X is 

the core wave number at a vacuum wavelength A. The evolution equation for the envelope 
is the paraxial wave (or paraxial Helmholtz) equation which reads 

2ikoU^ + Upp + p~^Up + p~‘^Up^ - V(p)u = 0 (2) 

where V(p) = —2klA{p/ay. The above equation is often viewed from a quantum- 
mechanical perspective namely as Schrodinger equation for wave function u{p, p, z) with ^ 
playing the role of time. In this context an antiguide is equivalent to an unstable poten¬ 
tial V (p) which justihes the minus sign that we used in its dehnition. We now assume a 
non-dihracting solution of Eq. ([2]) of the form 

u = U{p) exp {imp + ik^bz) (3) 


where integer m is the order of vorticity (or topological charge) and 6 is a dimensionless real 
number with |6| 1 so that the assumption of a slowly varying along z envelope stays true. 

Substituting into Eq. ([2]) we get 


p^U"{p)+pU'{p) + 


- 2hky - m^ ) U{p) = 0 


( 4 ) 


with the primes indicating ordinary differentiation with respect to p. The above falls under 
the general family of Bocher’s equations 23| and has no tabulated solution for a general 
value of V. The solution has to be obtained numerically or in the form of a series, a method 
that is frequently applied to graded-index optical hbers. [2^ In the particular case v = 2 and 
6 < 0, Eq. dl]) is called the Bessel wave equation and has been studied to some extent in 
the literature. In analogy to the terminology of Bessel functions, its non-singular solution is 
termed Bessel wave function of the hrst kind JJm { ko\/2A/a, ko\/—2b,py 25| 
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Let us first note two cases where Eq. (jl]) admits a solution in terms of standard Bessel 
functions. The hrst is the limit A = 0 whereby the medium is homogeneous with an index 
n = riQ. Then Eq. (jl]) reduces to Bessel equation with the nonsingular solution being the 


Bessel function of the hrst kind 


26| 


U (p) = Jm pj 


(5) 


We thus obtain a standard Bessel beam of order m. Obviously we must have 6 < 0, which 
means that the propagation constant kz = fco(l + &) is lower than the medium wave number 
fco (recall the factor muliplying the envelope function u of Eq. (j^). 

The second case is that of a wave with 6 = 0 namely with a propagation constant kz = ko 
equal to the core wave number. We then obtain a transformed version of Bessel equation 
with nonsingular solution 


U{p) 



( 6 ) 


This is a chirped Bessel beam with vorticity m that propagates without spreading (or 
diffracting) along the antiguide. Note that, for real u and integer m, the order of the 
Bessel function is generally a real number which reduces to m only when u = 0 namely 
in a homogenous medium. Why such a solution is chirped we can easily understand from 
the large-argument asymptotics of Bessel functions 


Jn{x) ~ 



(7) 


which holds generally for real order n (and complex argument with |arp(a;)| < vr). 

Having a closed-form solution like (|6]) at hand, few interesting remarks are in order. First, 
note that the argument of the Bessel function in Eq. ([6]) can be intuitively explained by the 
concept of the local radial wave number k{p). This is a measure of the spatial frequency of 
the radial oscillations and, due to the asymptotic form (I7|), it is dehned as the derivative of 
the argument in Eq. (|6]) with respect to p. The result is 

k (p) = p-'" = y v/ToTTg (8) 


where we have also used Eq. ([T]). Thus k{p) and kz are connected through the dispersion 
relationship of a plane wave in a homogeneous medium with index equal to the local value 
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n{p). Second, note that in a homogeneous medium with index n = no, ko is a cut-off 
propagation constant for propagating waves in the sense that kz < k^. As kz —)■ k^ , the 
transverse wave number tends to zero and the corresponding Bessel beam oscillates so slowly 
that it effectively assumes a constant amplitude which is nonzero only for the Jq case. 
However, when the index is perturbed as in Eq. ([T]), waves with kz > k^ are enabled 
and, no matter how small A is, chirped Bessel waves can oscillate many times as long as the 
characteristic radius is large enough compared to the wavelength in the medium (a S> A/no). 
This is easily deduced from Eq. ([6]) by setting p = a. The argument becomes proportional 
to fco®\/A thus allowing a large ratio ano/X to compensate for a small A. 


It is also interesting to see how the solution of Eq. (jH) behaves near the axis of the 
antiguide. In the limit p —)■ 0, the power term becomes negligibly small compared to 
the quadratic term and Eq. (0]) reduces to the Bessel or modihed Bessel equation depening 
on whether 6 < 0 or 6 > 0, respectively. Thus in this limit the solution behaves either like 
a Bessel or modihed Bessel function of the hrst kind with the same argument with Eq. ([5]). 
This could also be expected intuitively since, near the axis, the medium has a hat index 
prohle n ^ no hence the wave should behave as a standard Bessel beam in a homogeneous 
medium. In any case the small argument asymptotic of Bessel functions reads 261 


Jn {x) ~ 


T (n -h 1) 


(9) 


hence the envelope of the wave behaves near the axis as u ct. X In the other 

extreme, namely as p —>■ oo, the power term dominates the quadratic term which can 
be neglected. We thus obtain the modihed Bessel equation of the 6 = 0 case and, since we 
are far from p = 0, the solution behaves asymptotically as a linear superposition of a Bessel 
and a Neumann function (or Bessel function of the second kind) with the argument of Eq. 


Figured] shows examples of chirped Bessel waves with 6 = 0 and various orders of vorticity. 
For the antiguide we have used parameters that are common to graded-index optical hbers. 
One can observe the chirped oscillations of the solution as well as the discussed behaviour 
near the axis. In the m = 0 case notice the battening of the amplitude near the axis which 
is due to the Jo(ct. x p^) dependence. 
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FIG. 1. (a) Amplitude U{p) of the chirped Bessel waves of Eq. ([6]) for different orders of vorticity 
m = 0,1,2,3. The insets show the corresponding transverse intensity profiles \U\^. For the 
antiguide we have chosen the parameters no = 1.5, A = 0.01, u = 2, a = 50pm at a vacuum 
wavelength A = 1.5^m. 


B. Ray analysis 


Before proceeding to numerical simulations, it is instructive to obtain a better physical 
intuition of the waves under consideration using ray optics. We have already mentioned 
that the non-diffracting nature of the ID Airy wavepacket in a linear potential (constant 
force held) can be explained through its structure of parabolic rays enveloping a straight 
caustic. A generalization of this picture to a nonlinear power-law potential provides the ray 
structure of our chirped Bessel waves. For simplicity we limit ourselves to waves with zero 
angular momentum (m = 0). Then the rays are meridional and the paraxial ray equation is 
written in the p — z plane as 2^ 

cPp dn (p) 


n{p) 


dz"^ 


dp 


( 10 ) 


where we have used the fact that the index of refraction does not dependent on z. As a 
result, the ray equation is invariant along the axis or, using the strict mathematical term, 
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autonomous. This implies that, if p{z) is a solution to Eq. (1T0|) . the shifted along the z axis 
ray p{z — L) is also a solution for any L. Based on this property, the ray pattern of a non¬ 
diffracting wave is obtained as a continuum of shifted copies of a particular solution of the 
ray equation for L running ideally from — oo to cxo. By revolving this p — z pattern around 
the z axis, the total 3D ray structure is obtained. A transverse cut of this ray structure at 
any level ^ yields the same interference pattern which is what one physically observes as a 
non-diffracting wave. 

By substituting Eq. ([1]) into the ray equation flTOD and rearranging we get 


d^p _ Au 
dz‘^ a'' ^ 


( 11 ) 


where in the right-hand side we have approximated a factor [1 -|- 2A {p/aY]~^ by unity. For 
V = 1 {a. linear potential), the right-hand side is a constant and we obtain parabolic rays of 
the form 

Pi (z) = + ^2 (12) 

where C'i,C '2 are integration constants while the subscript I is a reminder of the linear 
variation of the refractive index. For u = 2 (inverse harmonic oscillator potential), the ray 
equation is solved by exponential rays of the general form 


Pq {z) = Cl exp 



-h C 2 exp 



(13) 


where the subscript q indicates the quadratic variation of the refractive index. For z/ 7 ^ 2, 
Eq. fllip is a second-order nonlinear differential equation of the Emden-Fowler type 23| and 


has a known implicit solution. 


28| Instead of directly applying the relevant formula, we 


here hnd it more useful to reduce Eq. flTT]) to an equivalent hrst-order differential equation. 
This is done by multiplying both sides with dp/dz and integrating with respect to z. The 
result is 


f, - 

where C 3 is a dimensionless constant of integration. This constant is connected with the 
propagation constant of our non-diffracting waves. This is better understood if one views 
Eq. flTTl) as the conservation law of the longitudinal momentum in our system. Indeed, in a 
medium that is invariant along z, the component of the wavevector of a ray is conserved 
exactly as the tangential wave vector after refraction through a planar interface (Snell’s law). 
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It follows that each ray is characterized by a kz that stays invariant along its trajectory and 
is given by the projection of the local wave vector along the 2 ; axis, namely 


kz 


2ti 

T 


n (p) cos 9z 


271 


n{p) 



-1/2 


(15) 


where dz = arctan (dp/d^) is the angle between the ray and the z axis. According to our 
dehnition (|3]) and the previous discussion, the rays that compose our non-diffracting waves 
have the invariant kz = ko{l + b). Substituting into Eq. flT^ and making a paraxial 
approximation (|dp/d^| <C 1), we get 

Tz = (“> 

which is identical with Eq. flT^ if = —2b. 

The particular form of the ray equation ffTTl) or flTbl) allows to distinguish two different 
regimes or patterns of ray structure. For & > 0, the slope of a ray with respect to the axis 
dp/dz vanishes at the radius Pc = a ■ This radius is a turning point of the rays and 

is identihed with a caustic radius of the corresponding non-diffracting wave. Since the rays 
cannot reach this region, the intensity of the wave is low in p < pc- For waves with b < 0, 
on the other hand, dp/dz never vanishes and the rays cross the axis of the antiguide with 
a slope ±\/—25. After the crossing their convexity changes sign and the total trajectory is 
symmetric with respect to the crossing point. In the boundary between the two regimes, 

6 = 0 and the ray equation can be integrated analytically to yield 

2 

2-z/ 

(17) 


Pfe=o {z) = a 


( 


V 


\2 


-1 (L±z) 


where the sign is chosen to yield a positive number inside the brackets. Thus the rays either 
touch the z axis with zero slope when u < 2, either they are asymptotic to the axis when 
V > 2. In the latter case the rays are also asymptotic to transverse planes z = ±L. Finally, 
for 6 = 0 and u = 2, Eq. flTT)) yields exponential rays which are also asymptotic to the 2 ; 
axis. These rays are also obtained from Eq. flT^ for Ci = 0 or (72 = 0. 

Figure [2] illustrates the different cases of rays discussed above. The rays are repeated 
along the 2 ; axis to show the ray structure of our chirped Bessel waves in a meridional plane. 
The slope of the rays has been exagerrated for visualization purposes. Figure [2](a) shows 
the case 6 > 0 where the rays are reflected at the caustic radius Pc which is the point at 
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FIG. 2. Ray patterns of non-diffracting waves in an optical antiguide for (a) b > 0, (h) b < 0, (c) 
6 = 0 and i' < 2, (d) 6 = 0 and v >2. The red curves on the right of the patterns are the relative 
index difference 5n{p) = (n(p) — no)/no with the value of 6 indicated. In (b) and (d) two different 
ray colors are used for visualization purposes. 


which b equals the relative index difference bn{p) ~ A{p/ay. In (b) we have the case 6 < 0 
where the rays cross the axis and continue with the opposite convexity. In (c) we have the 
case 6 = 0 and v < 2 m which the rays are reflected at the axis, while (d) shows the case 
6 = 0 and v > 2 where the rays are asymptotic to the axis. Finally note that these ray 
patterns suggest that the chirped Bessel waves are also able to self-heal like their standard 
Bessel counterparts, since the rays contributing to the held near the axis at a large distance 
z are emitted from far off-axis points on the input plane z = 0 . 

The above concern waves without angular momentum. By introducing the azimuthal 
coordinate in the ray equation, ray patterns of vortex waves can also be derived. These 
involve skew rather than meridional rays. The vorticity introduces an additional unstable 
potential with a dependence w?/p^ which acts as a centrifugal force pushing the intensity 
distribution of the waves away from the axis which becomes a held node. Finally note that, 
in using paraxial ray optics results, one should be aware of their limitations. A rule of thumb 
is that the paraxial approximation is satisfactory for ray angles smaller than 30° or roughly 


slopes \dp/dz\ < 0.5. For larger angles, the cornp 


to determine the accurate shapes of the rays. 


ete nonparaxial ray equation is required 
24| For example, rays in a medium with a 
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linear index profile are actually parabolic only within a finite propagation length around 
their turning point, where the paraxial approximation holds, while they diverge with an 
exponential law at larger slopes. However, we will not pursue ray optics any further here. 


III. NUMERICAL EXPERIMENTS 


We now proceed to investigate numerically the propagation of chirped Bessel waves in 
antiguides. As mentioned, their radial prohle is obtained by solving Eq. (jl]). To this 
end, we transform this second-order differential equation to a system of two first-order 
differential equations for functions U{p) and U'{p) and then we solve this system using a 
high-order Runge-Kutta scheme. 291] To avoid the singularity of this system at p = 0, the 
initial conditions are imposed at p = e > 0 and their values are obtained from the limiting 
behaviour of the solution that we determined in section III Al For a wave with m = 0 for 
example, we may start with U{e) = 1, U'{e) = 0. 

Figures [3]]^a) and (b) show the radial prohle of chirped Bessel waves with m = 0 and 
6 < 0 (^2 < fco) or 6 > 0 {kz > ko), respectively. The parameters of the antiguide and the 
operating wavelength are the same with Fig. [H For 6 < 0 the wave resembles a zero-order 
Bessel beam with a maximum on axis and oscillations that are getting faster with increasing 
p. For 6 > 0, the maximum shifts away from the axis which rehects the existence of a 
caustic radius, namely the turning point of the rays that was discussed in Fig. In 

both cases the oscillations are linearly chirped (phase increasing with p^) as expected from 
Eq. dH]). For comparison, the Bessel functions that are asymptotic to the solutions as p —)■ 0 
have been added to both figures. These are given by Eq. (jS]). When 6 > 0 the latter 
equation suggests a modihed Bessel function. The complete evolution of the solution with 
b is tracked in the map of Fig. He). As b assumes lower negative values, the oscillations 
get compressed because the local transverse wave number increases (recall the discussion 
following Eq. ([8]))- As 6 assumes larger positive values, the maximum shifts further away 
from the axis because the caustic radius increases. Exactly as happens with the Airy wave 
packet, the wave amplitude decays super-exponentially toward the dark side of the caustic 
and thus the interior of the caustic cylinder contains low optical intentisy. 

We can now use the chirped Bessel waves that we have determined numerically as an input 
to an optical antiguide and investigate their propagation dynamics. At this point it should 
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FIG. 3. Radial profiles (blue curves) of chirped Bessel beams obtained by solving Eq. ([1]) numeri¬ 
cally with (a) b = —0.001 {kz = 0.999A:o) and (b) b = 0.001 {kz = l.OOlfco). The thin red curves 
are the Bessel functions that are asymptotic to the solution as p —)• 0 given by Eq. m- (c) A 
continuous map of the magnitude |17| versus p and b. For any b the maximum intensity has been 
scaled to unity. The dashed white curve is the caustic radius pc = a(6/A)^A, antiguide’s 

parameters and the operating wavelength are the same with Fig. [1] 

be realized that, similar to standard Bessel beams in homogeneous media, the solutions 
of Eq. dH) carry inhnite power and thus cannot be realized in their ideal mathematical 
form. To see this, consider the integral of U^{p)n{p)pdp which is proportional to the power 
carried by the wave. Then, as discussed in the asymptotics of Section Hi Al for increasing p, 
U{p) ~ p-(^+2)/4 -which makes this integral diverge. Hence one needs to consider realistic, 
hnite-energy counterparts which are obtained by a smooth or abrupt apodization of the ideal 
solutions. 
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FIG. 4. (a) Propagation dynamics in terms of the intensity \U{x,0, z)\‘^ on the y = 0 plane for a 
chirped Bessel wave with m = 0 and b = —0.001 (radial profile shown in Fig. EKa)), apodized with 
the Gaussian exp[—(/9/125)^] (all distances in /rm), in an antiguide with the parameters of Fig. 
[TJ (b) Amplitude profile \U{x,0,z)\ of the wave of (a) at different 2 ; levels, (c-d) Corresponding 
results for the Gaussian input wave U{p,z = 0) = exp[—(/?/4)^]. For visualization purposes, (c) 
shows the amplitude \U\ instead of the intensity \U\‘^. 


Figure 111(a) depicts the propagation of a chirped Bessel wave with zero angular momen¬ 
tum [m = 0) and a given kz in a quadratic antiguide. The input condition is apodized with 
a Gaussian function. The simulation is a split-step Fourier numerical integration of Eq. (E]). 


30l | The result clearly verihes the non-diffracting character of the wave, which is manifested 


in its apodized version as a remarkable resistance to diffraction. This becomes evident in 
more detail in Fig. l4](b) which shows the radial distribution of intensity at various distances 
along the antiguide. Notice how the prohle, and particularly the central lobe, persists with 
increasing z despite the defocusing action of the unstable refractive index prohle. As visu¬ 
alized with ray optics in Section III Bl the outer intensity rings constantly “feed” the central 
lobe with energy thus compensating for the spreading effect of the antiguide. Moreover, the 
interference of these contributions is such that its shape remains unchanged over an appre¬ 
ciable propagation distance which approaches 1 mm for the given apodization function. To 
appreciate the uniqueness of this effect even more, in Figs. 111(c) and (d) we consider the 
propagation of a Gaussian input wave with the same full width at half maximum (FWHM). 
This is similar to launching just the central lobe of the chirped Bessel wave. In the absence of 
some special structure in its wavefront, such a wave has no chances to survive the defocusing 
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FIG. 5. (a) Propagation dynamics in terms of the intensity \U{x,0, z)\‘^ on the y = 0 plane for a 
vortex chirped Bessel wave with m = 1 and b = —0.001, apodized with the Gaussian exp[—(/9/125)^] 
(all distances in fim), in an antiguide with the parameters of Fig. [H (b) Amplitude profile 
\U{x,0,z)\ of the wave of (a) at different z levels, (c-d) Corresponding results for the Laguerre- 
Gaussian input wave u{p,(p,z = 0) = pex.p[iip — (p/7.5)^]. For visualization purposes, (c) shows 
the amplitude \U\ instead of the intensity |C/p. 


action of the antiguide and thus spreads quickly. These results already suggest that chirped 
Bessel waves can outperform standard waves (like Gaussian or the simple Laguerre-Gaussian 


modes used in 


2l| ) in applications where light has to be guided in a diffraction-resisting 


manner in unstable optical potentials. 


Figure [5] shows another example, this time for a vortex wave with m = 1. Subhgures (a) 
and (b) demonstrate the resistance of this wave to diffraction in the same antiguide with 
Fig. [U By contrast, in (c) and (d) we observe the strong spreading of a wave with input 
amplitude of the form pexp{iip — namely the focal-plane expression of the lowest 

Laguerre-Gauss mode with the same order of vorticity. Again, the structure of the chirped 
intensity rings is what keeps the chirped Bessel wave undistorted along the antiguide, in 
contrast to the single ring of the Lagueree-Gauss mode which spreads out quickly. 


Let us now test our waves under some abnormal conditions. The hrst such condition is 
eccentric launching and is examined in Fig. Ej where the chirped Bessel wave of Fig. IHis 
launched with its axis displaced by few pm from the axis of the antigide. The result is quite 
remarkable: The wave is attracted toward the high-index region along a curved trajectory 
in the direction of the displacement [x in this example) while retaining its transverse profile. 
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FIG. 6. (a) Propagation dynamics in terms of the intensity \U{x,Q, z)\^ on the y = 0 plane for 
the wave of Fig. 0] launched with a 8 ym offset from the antiguide’s axis in the x direction, (b-c) 
Transverse intensity profiles at the shown values of z. (d) Amplitude profile \U{x, 0, z)\ of the wave 
of (a) at different z levels. 


In the amplitude profiles shown in (d) one can indeed verify the impressive robustness of 
the main lobe during propagation. 


The second adverse condition is launching the wave in a curved antiguide. To emulate 
the bending effect we solve Eq. ([2]) with a 2 ;-dependent potential of the form Vb{x, y, z) = 
V where V is the potential dehned below Eq. ([2]) and f(z) is the 

function describing the curve of the antiguide, the simplest choice being a parabola, i.e. 
f{z) = 'jz'^. Chirped Bessel waves of zero (a-b) and nonzero (c-d) angular momentum are 
examined. In both cases the wave starts propagating along a straight trajectory demonstrat¬ 
ing its non-diffracting property. After some distance (around 400 fim in both examples) and 
due to the bending of the antiguide, the wave hnds a large part of its wavefront displaced 
with respect to the antiguide’s axis. Thus, as in the case of eccentric launching, its tra¬ 
jectory bends toward the region of increasing index (to the left in these examples). Still, 
despite this deformation of trajectory the amplitude prohle remains fairly robust. Notice 
the difference of these results with the similar numerical experiment reported in 2l|] using a 
simple Laguerre-Gauss mode with vorticity m = 1. In the latter case and due to the absence 
of some special structure in its wavefront, the wave gets severely distorted. 
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FIG. 7. (a) Propagation dynamics in terms of the intensity \U{x,Q, z)\‘^ on the y = 0 plane for the 
wave of Fig. IHlaunched in an antiguide whose axis is curved along to the parabola f{z) = (z/133)^ 
(white dashed line, all distances in /rm). (b) Amplitude profile \U{x,^,z)\ of the wave of (a) at 
different 2 : levels, (c-d) Corresponding information for the vortex wave of Fig. [5] launched in the 
same curved antiguide. 


The third and last adverse scenario aims to demonstrate self-healing which we have 
already mentioned as a property of the chirped Bessel waves while discnssing their ray 
optics. An example is shown in Fig. [H] where the chirped Bessel wave of Fig. IHis lannched 
in the same qnadratic antignide with its main lobe removed. Despite this severe distortion, 
the wave manages to reconstruct the main lobe and part of the profile after some distance. 

We close with a more practical approach to the above effects. In practice, it may be 
difficult to realize an exact wave function experimentally. This might be true not only for 
our chirped Bessel waves but for any kind of structured lightwave. It is thus always useful 
to come up with simpler input waves able to approximate the same effect which in our case 
is the diffraction-resisting propagation in an optical antiguide. From a closer look at the 
ray patterns in Figs. [2l one sees that, at any given 2 ;, there are two rays passing from a 
point on the p axis, an inward and an outward one, with equal but opposite slopes with 
respect to the z axis. Of the two, it is the inward ray that contributes to the held near the 
axis at larger values of This implies that the prohle of the wave near the axis will not 
be affected signihcantly if we discard the outward rays from the input condition ai z = 0. 
Such a condition can be easily implemented in the laboratory with a circularly symmetric 
phase mask (or a spatial light modulator) that has the appropriate phase function <I>(p) to 
create the inward rays with the correct slope. *h(p) is obviously the phase of the oscillations 
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FIG. 8. (a) Propagation dynamics in terms of the intensity \U{x,Q, z)\‘^ on the y = 0 plane for the 
chirped Bessel wave of Fig. |3]with its central lobed blocked, (b) Amplitude profile \U{x,0,z)\ at 
different z levels for the wave of (a) (red lines) and the wave of Fig. H^a) (blue lines). 

of the solution of Eq. (01) (assuming that it oscillates proportionally to cos(<h(p)) and the 
transmission function of the phase mask should be exp[—i$(p)], with the negative sign 
giving rise to inward rays. An relevant example is given in Fig. |9l where we are trying to 
approximate the dynamics of Fig. 0](a) with a Gaussian input wave modulated by the phase 
mask shown in (c-d). For a better comparison, (b) shows the amplitude prohle at several 
values of the propagation distance superposed on the corresponding prohles of Fig. 111(b). 
The goal is obviously achieved as the phase-modulated Gaussian wave reproduces almost 
excellently the main lobe and the first intensity ring of the chirped Bessel wave, within a 
certain range of z. 


IV. CONCLUSION 

We have introduced and studied chirped Bessel waves capable of resisting diffraction in 
optical antiguides with a power-law varying index of refraction. We examined the ray struc¬ 
ture of these waves with a detailed (almost tutorial) spirit and investigated their propagation 
dynamics with full-wave split-step Fourier simulations. Conditions of eccentric illumination. 
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FIG. 9. (a) Propagation dynamics in terms of the intensity |[/(x,0,z)p on the y = 0 plane for 
the input wave exp[i<I>(/9) — (p/60)^] launched in the antiguide of Fig. [D (all distances in ym). (b) 
Amplitude profile \U{x,0,z)\ at different z levels for the wave of (a) (red lines) and the wave of 
Fig. UK a) (blue lines), (c) The modulo-27r phase function ^{x,y) that modulates the Gaussian 
wave, (d) The same function (in radians) versus the radius. 


antiguide bending and partial obstruction of the input wavefront were also examined and 
demonstrated the robustness of their prohle and the ability for self-healing. A practical 
implementation of these waves using simple phase-modulated input waves was hnally pro¬ 
posed. From a more general viewpoint, this work adds to the family of structured lightwaves 
members that resist (or completely overcome, in their ideal mathematical version) diffrac¬ 
tion in media that are inherently more diffractive than free space, such as unstable optical 
ducts. And we have done so by extending the family of Bessel beams with their radially 
chirped counterparts. The former are non-diffracting in free space while the latter in optical 
antiguides. 
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